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ABSTRACT 

The active galaxy NGC1275 lies at the centre of the Perseus cluster of galaxies, which is 
the X-ray brightest cluster in the Sky. The nucleus shows large variability over the past few 
decades. We compile a lightcurve of its X-ray emission covering about 40 years and show 
that the bright phase around 1980 explains why the inner X-ray bubbles were not seen in the 
images taken with the Einstein Observatory. The flux had dropped considerably by 1992 when 
images with the ROSAT HRI led to their discovery. The nucleus is showing a slow X-ray rise 
since the first Chandra images in 2000. If it brightens back to the pre-1990 level, then X-ray 
absorption spectroscopy by ASTRO-H can reveal the velocity structure of the shocked gas 
surrounding the inner bubbles. 
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1 INTRODUCTION 

The Perseus cluster of galaxies, the X-ray brightest cluster in the 
Sky, has an X-ray cool core surrounding the active central galaxy 
NGC 1275. The jets forming its radio source 3C84 blow 15 kpc 
diameter bubbles in the hot gas (Boehringer et al 1993; Fabian et 
al 2000; 2006). Buoyant ghost bubbles lie further from the nucleus, 
revealing the last 10 8 yr or more of energy feedback from the AGN 
to the cool core. 

Radio and 7 —ray light curves of the active nucleus of 
NGC 1275 show large amplitude variability over the past 40 years 
(Dutson et al 2014; Abdo et al 2009). The source was bright from 
the 1960s to about 1990 when it abruptly faded by about an order 
of magnitude. It was faintest around 2000 and now shows signs of 
increasing again. 

Here we look at the long-term variability of NGC 1275, finding 
it to have dropped by about a factor or 20 between when first im¬ 
aged with the Einstein Observatory in 1980 and Chandra in 2000, 
rising slowly since then. We demonstrate that the high brightness 
of the nucleus prevented detection of the inner bubbles with the 
Einstein Observatory in the (relatively) short exposure times used. 

We also examine the possibility of detecting resonance line 
absorption from the hot gas using the nucleus as a backlight. It ap¬ 
pears that if the the nucleus remains at the 2006 level or similar 
then such absorption would be challenging with ASTRO-H, due 
to its half energy width of 1.5 arcmin causing the nucleus to ap¬ 
pear at low contrast against the bright hot gas emission. It will be 
straightforward for Athena with its angular resolution of 5 arcsec 
or better. If, however, the nucleus increases to the 1980 level, then 
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it will compromise emission line studies of the central region with 
ASTRO-H, but make absorption spectroscopy possible. Since our 
line of sight passes directly through the shocked gas surrounding 
the 1100 km s -1 expanding Southern bubble, such studies would 
give unique insight and test the growth and energetics of the bub¬ 
bles. 


2 OBSERVATIONS 

We investigate the variation in X-ray flux from the central AGN 
in Perseus by analysing all of the spatially resolved X-ray obser¬ 
vations that have been taken of the system, starting with the first 
Einstein observation from 1979. The observations used are listed 
in Table Q] In addition we have compiled a comprehensive list of 
hard energy spectral measurements for NGC 1275 from Coperni¬ 
cus, OSO-7, HEAO-1 A-4, EXOSAT LE, Ginga LAC, OSSE and 
several balloon experiments (Table 1). FWHM for the imaging in¬ 
struments is tabulated. The early Copernicus estimate contains an 
uncertain level of cluster emission so is shown as an upper limit. 

From the cleaned events files we find the count rate within a 
circular region centred on the AGN of diameter equal to the Full 
Width at Half Maximum (FWHM) for the instrument used. The 
count rate was then converted into the 0.5-10.0 keV flux using 
WebPimms, under the assumption that the spectrum of the AGN 
is a powerlaw of index 1.65, as was found in the XMM-Newton 
study of Churazov et al. (2003). Residual cluster emission is un¬ 
likely to contribute significantly to the measured nucleus flux when 
it is bright. 

The resulting X-ray lightcurve is shown in Fig. [TJ with points 
from different missions labelled and shown with different colours. 
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Figure 1. Variation in X-ray flux from the central AGN in Perseus from the 1970s 
by Dutson et al. (2014), which follows a similar trend to the X-ray lightcurve. 


The X-ray lightcurve shows a decrease by nearly an order of mag¬ 
nitude from the early 1980s to the 1990s. This is in good agreement 
with the decrease observed in the 90GHz radio emission reported 
in Dutson et al. (2014), which is plotted as the grey points. 


3 SIMULATED EINSTEIN IMAGES 

The Einstein Observatory’s HRI instrument observed the Perseus 
cluster for 15.6ks in 1979 February, and again for 6.6ks in 1980 
February. The central cavities cannot be resolved in either of these 
images, even when they are stacked togetheiQ. The cavities re¬ 
mained undiscovered until the first ROSAT HRI observations of 
Perseus (Boehringer et al. 1993), which showed for the first time 
clear evidence for the interaction between the relativistic particles 
in the radio lobes and the intracluster medium (Fig. 2). 

The on axis spatial resolution of the Einstein HRI was, how¬ 
ever, comparable to that of the ROSAT HRI (3 arcsec FWHM for 
Einstein, 2 arcsec FWHM for ROSAT). It is therefore interesting to 
see whether the high X-ray flux from the central AGN during the 

1 A decrease in X-ray emission at the position of the the NW ghost cavity is 
apparent in the Einstein image, and was commented on at the time (Fabian 
et al 1981, Branduardi-Raymont et al 1981). Since it is devoid of detected 
radio emission, its true nature was unsuspected. 


to the present day. The grey points show the 90GHz radio lightcurve gathered 


Einstein mission prevented Einstein from making the first detection 
of the cavities in Perseus. 

We therefore simulated Einstein images of Perseus to investi¬ 
gate the visibility of the central cavities as a function of the central 
AGN flux and exposure time. The simulations were performed us¬ 
ing SIM30, which was modified to allow it to simulate the Einstein 
observatory. This was achieved by setting the focal length, field 
of view and pixel size parameters equal to those for the Einstein 
HRI, and using Einstein RMF and ARF response files. The PSF 
was set by using the Encircled Energy Function from Van Spey- 
broeck (1979). For the input image we used the stacked, 890ks, 
Chandra image from Fabian et al. (2006). Since only one spectrum 
can be used when simulating an events file with SIMX, we simulated 
separate events files for the central AGN and the cluster emission 
with different spectra for each, and then added them. For the central 
AGN the input spectrum was set as a powerlaw of index 1.65, as 
found in Churazov et al. (2003). For the cluster emission the input 
spectrum was a thermal APEC component with parameters fixed to 
the mean cluster temperature in the core (kT =4keV). 

We perform Einstein HRI simulations with 3 different AGN 
flux levels and 3 different exposure times, and these are shown in 
Fig. El The top row in Fig. [3] shows simulations using the 1979 

2 http://hea-www.cfa.harvard.edu/simx/ 
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Figure 2. Clockwise from top left: Einstein HRI image of NGC1275; ROSAT HRI image of NGC1275; Einstein HRI image of the bright quasar 3C273 to 
same angular scale; Chandra image of NGC1275. The Einstein image has been lighly smoothed. 


AGN flux level (7.3 x 10 -11 erg s -1 cm -2 in the 0.5-10keV band) 
for exposure times of lx, 10x and 100x that of the 15.6ks Ein¬ 
stein observation from 1979. With the AGN at its 1979 flux level, 
the PSF spreading prevents the central cavities from being robustly 
measured even in very deep (1.5Ms) observations. The north west¬ 
ern ‘ghost’ cavity is however resolved in the ^156ks simulations. 

The second row shows the case when the central AGN flux 
is reduced to its minimum observed level (from around 2002) of 
5xl0 -12 erg s -1 cm -2 in the 0.5-10 keV band. We see that a 
15.6ks observation can begin to resolve the inner cavities (1st col¬ 
umn, 2nd row in Fig. [3j. With ten times this exposure time, the 
cavity system is clearly resolved, showing a similar level of detail 
as the ROSAT HRI observations. 

Finally, in the third row, we show a simulated Einstein im¬ 
age of the case where there is no contribution from a central AGN 
(bottom row), allowing us to demonstrate Einstein’s ability to re¬ 
solve fluctuations in extended X-ray emission from the ICM in the 
idealised case of no point source contamination. We see that even 
in a 15.6ks observation, matching the exposure time of the 1979 
Einstein observation, the inner cavities can just about be resolved. 
Structure is seen on the scale of the radio source which may at least 


have triggered further observation. The cavity system is comfort¬ 
ably resolved with an exposure time which is ten times deeper. 


4 X-RAY ABSORPTION SPECTROSCOPY 
USING NGC1275 

In this section, we analyse absorption by the ICM acting on the 
AGN emission, using a multiphase collisionally-ionized absorber. 
To calculate the column density of gas in each temperature range, 
we extracted azimuthally averaged temperature and density profiles 
from the Chandra data, which were deprojected using the DSDE- 
PROJ direct spectral deprojection method (Sanders & Fabian 2007; 
Russell, Sanders & Fabian 2008). The reduction of the Chandra 
data is described in Sanders & Fabian (2007). We integrated the de¬ 
projected density profile along the line of sight to find the column 
density in the following temperature bands which are tabulated in 
Table [5] <2, 2-3, 3-4, 4-5 and 5-6 keV. 

The lowest two temperature components are associated with 
the high pressure (shocked) rims around the bubbles (Fig. 6). The 
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Table 1. Archival Perseus cluster observations 


Telescope 

Detector 

FWHM Obs ID 

arcsec 

Year 

Day 

Exposure Reference 
ks 

Copernicus 


120 


1972 Sep 



Fabian et al. (1974) 

UCSD OSO-7 




1972 9-27 Nov 



Rothschild et al (1981) 

HEAO-1 A4 




1977 Aug; Feb, 1978 Aug 



Primini et al. (1981) 

Einstein 

HRI 

3 

h0316n41.xia 

1979 

52 

15.6 

Fabian et al. (1981) 

Einstein 

HRI 


h0316n41.xib 

1980 

49 

6.6 

Fabian et al. (1981) 

TESRE Balloon 



1981 29 Sep 



Matt et al. (1990) 

EXOSAT 

LELX3 

15 

ex830724 

1983 

205 

13.4 

Allen et al. (1992) 

EXOSAT 

LE LX3 


ex840124 

1984 

24 

2.3 

Allen et al. (1992) 

SPARTAN-1 




1984 Aug 



Ulmer et al. (1987) 

SPACELAB-2 




1985 Sep 



Eyles et al. (1991) 

Ginga 




1989 

15-17 


Allen et al. (1992) 

OSSE 




1991 28 Nov-12 Dec 



Osako et al. (1994) 

ROSAT 

HRI 

2 

rh800068n00 

1991 

41 

10.9 

Boehringer et al. (1993) 

ROSAT 

HRI 


rh800591a01 

1994 

217 

52.9 


ROSAT 

HRI 


rh702626n00 

1996 

42 

8.2 


ROSAT 

HRI 


rh702626a01 

1997 

62 

4.6 


XMM 

PN 

7 

P0085110101PN 

2001 

30 

51.2 

Churazov et al. (2003) 

XMM 

PN 


P0305780101PN 

2006 

29 

119 


Swift 

XRT 

7 

sw00036524002xpcw2po 

2007 

340 

3.6 


Swift 

XRT 


sw00030354003xpcw3po 

2009 

364 

4.3 


Swift 

XRT 


sw00031770009xpcw3po 

2010 

219 

2.1 


Swift 

XRT 


sw00049799006xpcw3po 

2013 

213 

1.6 



Table 2. Copernicus: 12, 6, and 2 arcmin apertures; OSO-7: 6.5 deg FWHM; HEAO-1: 1.7 x 20 deg; TESRE balloon: 3 deg; SPARTAN-1: 5 arcmin x 3 deg; 
SPACELAB-2: 12 arcmin; Ginga: 0.8 x 1.7 deg; OSSE: 3.8 x 11.4 deg. 
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Figure 3. Simulated Einstein HRI images with different levels of central AGN X-ray flux and different exposure times (exposure times are 1 x, 10x and lOOx 
the 1979 Einstein observation exposure of 15.6ks). 
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Table 3. Multiphase absorption model 


Trange (keV) 

Tadopted (keV) 

A/H(10 21 cm 2 ) 

< 2.0 

1.0 

2.3 

2.0-3.0 

2.5 

5.0 

3.0-4.0 

3.5 

4.0 

4.0-5.0 

4.5 

2.8 

5.0-6.0 

5.5 

1.0 


Notes: the model consists of five isothermal ionized absorber in collisional 
equilibrium, each of which has fixed temperature T ac j op ted an d hydrogen 
column density Nu . 



Figure 4. Transmission of the multiphase absorption model. 

density jump at the outer edge of the rims is 1.3 and the inferred 
Mach number is 1.2 (Graham et al 2008). 

For the absorption, each phase was modeled with a single¬ 
temperature hot model in SPEX0 (version 2.03.03), which cal¬ 
culates the transmission of a collisionally-ionized equilibrium 
plasma. For a given temperature and set of abundances, the model 
calculates the ionization balance and then determines all the ionic 
column densities by scaling to the prescribed total hydrogen col¬ 
umn density. We used five hot components with temperatures and 
hydrogen column densities as quoted in Table[3] We adopted a 
metallicity of 0.65 (Sanders & Fabian 2007) and a nominal ve¬ 
locity dispersion of 300 km s -1 We show the transmission of this 
multi-temperature absorber near the Fe K complex in Fig.[4]with 
the strongest absorption line labeled. An AGN spectrum provided 
with both high spatial and spectral resolution would reveal a strong 
absorption feature at 6.58 keV. 

In order to estimate the nuclear emission of NGC1275, we 
have analyzed the longest on-axis exposure of the Perseus cluster 
taken with XMM-Newton/EPIC in 2006 (ID=0305780101). Data 
were reduced and corrected by solar flares with the SAS vl4 fol¬ 
lowing the standard XMM-SAS procedure. We extracted the nu¬ 
clear spectrum and subtracted the cluster spectrum following the 
procedure of Churazov et al. (2003). The cluster-subtracted nucleus 
spectrum was then fitted with a power law and a gaussian to model 
the nuclear Fe K emission line, both redshifted and absorbed by 
the Galactic foreground. EPIC spectra do not have spectral reso- 


3 www.sron.nl/spex 



Figure 5. ATHENA/X-IFU simulation for the inner 14” . 


lution high enough to resolve any narrow absorption lines of the 
foreground ICM. 

The X-ray integral field unit (X-IFU) aboard ATHENA is 
planned to have spatial (5 arcsec) and spectral (2.5eV) resolution 
at 6keV high enough to resolve the Fe K lines absorbed by the 
nucleus. We used the nuclear spectral model fitted with the EPIC 
spectra as a template model to simulate a lOOks X-IFU spectrum. 
On top of the AGN emission the multiphase absorber was tested 
with three different values of velocity dispersion v a (see Fig.0. 
The high-quality X-IFU spectrum will easily resolve the lines. 

The ASTRO-H Soft X-ray Spectrometer (SXS) with its 5 eV 
spectral resolution could potentially resolve absorption in the nu¬ 
cleus emission of NGC1275, but due to the large 1.5 arcmin HEW, 
the SXS central chip will include both the central AGN and the 
cluster emission from a region of about 1.5 arcmin width. In or¬ 
der to estimate the cluster emission, we extracted a further EPIC 
spectrum in an annulus of 0.75 arcmin outer radius and inner ra¬ 
dius equal to that of the nuclear region (14 arcsec, see Churazov et 
al. 2003). We modelled the cluster spectrum with a 4keV isother¬ 
mal collisionally-ionized emission model (CIE model in SPEX). 
The abundances of all the atomic species from carbon to iron were 
fixed to 0.65 similarly to the multiphase absorber. 

We use these EPIC spectral fits to finally define a hybrid model 
which includes the nuclear emission (power law and Fe K line) in¬ 
trinsically absorbed by the multiphase ICM plus the unabsorbed 
1.5 arcmin cluster emission (CIE), both corrected by redshift and 
Galactic foreground. We simulate an ASTRO-H/SXS lOOks expo¬ 
sure for two models with and without AGN intrinsic absorption. 
We could not distinguish between the models because the cluster 
accounts for more than 95% of the emission within 1.5 arcmin. 
A very long exposure above 500 ks could potentially reveal some 
weak evidence of absorption. However, as previously mentioned, 
the AGN is brightening in recent years and at some point it may be 
as bright as in 1980, which was 10 time brighter than in 2006. We 
have therefore simulated a 500ks ASTRO-H/SXS spectrum with 
three models of AGN absorption calculated for a combination of 
different velocity dispersion and line-of-sight velocity of the mul¬ 
tiphase absorber (see Fig.lTJ). The lowest (blue) line in this Figure 
includes a velocity offset of —300 km s -1 (the bubble is likely ris¬ 
ing buoyantly in the cluster potential) with a velocity dispersion 
of 300 km s -1 as a rough approximation of the expected velocity 
structure behind the weak shock along our line of sight (Fig. 6). 
The effects of absorption of the AGN are clearly noticeable. 

Resonance line scattering of the emission from the intracluster 
medium in the cluster is discussed in detail by Zhuravleva et al 
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Figure 6. Left: Pressure map of the region around NGC1275. The bubbles are seen to be surrounded by spherical high pressure rings of shocked gas (see e.g. 
Graham et al 2007). Right: schematic diagram showing how the bubbles probably lie along our line of sight, with the nucleus being viewed through the top 
part of the S bubble. 



Energy (keV) 


Figure 7. ASTRO-H/SXS simulation of the central 1.5’ diameter cluster 
emission on top of an AGN spectrum lOx brighter than in 2006. Three mod¬ 
els of AGN instrinsic absorption are shown for a combination of different 
velocity dispersion and line-of-sight velocity. 

(2013) and is not included in our modelling. We note that the largest 
red and blueshifted absorption components seen in the spectrum are 
likely to be due to the expanding inner bubbles. 

We note that the shape of the 6.4 keV neural iron line (shown 
in Fig. 5) may hold interesting information (Churazov et al 1998) 
about the velocity structure of the innermost massive molecular 
clouds surrounding the nucleus (Salome et al 2006; see also Jaffe 
1990 for HI absorption against the nucleus). 


5 DISCUSSION 

We have shown that the X-ray bright nucleus can play both a neg¬ 
ative and positive role in the detectability of features in the intra¬ 
cluster medium. When the nucleus of NGC1275 was very bright 


back in 1980 it prevented detection of the inner cavities. Our un¬ 
derstanding of AGN feedback in clusters would have begun much 
earlier had the nucleus been weak at that time. If the nucleus bright¬ 
ens to a similar level in the next decade then it will enable unique 
absorption spectroscopy of the high pressure expanding rim around 
the inner Southern cavity to be carried out. 

Large variability of the jetted nucleus emission is common in 
cool core clusters. HST-1, part of the inner jet of M87, brightened 
and faded in X-rays by over an order of magnitude between 2004 
and 2008 (Harris et al 2009). Russell et al (2013) find significant 
X-ray variability in the nucleus of A2052 and Hydra-A; Hogan et 
al (2015) find high-frequency radio variability in 4 out of 23 central 
cluster sources. As shown above, an X-ray bright nucleus can en¬ 
able X-ray absorption spectroscopy. The quasar HI821+643 at the 
centre of an X-ray luminous cluster at redshift z — 0.297 ( Fang et 
al 2002; Russell et al 2010; Walker et al 2014) is a potentially good 
target. 
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